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ABSTRACT 

Recent analytical and numerical work argue that successful relativistic Fermi acceleration re- 
quires a weak magnetization of the unshocked plasma, all the more so at high Lorentz factors. 
The present paper tests this conclusion by computing the afterglow of a gamma-ray burst 
outflow propagating in a magnetized stellar wind using "ab initio" principles regarding the 
microphysics of relativistic Fermi acceleration. It is shown that in magnetized environments, 
one expects a drop-out in the X-ray band on sub-day scales as the synchrotron emission of 
the shock heated electrons exits the frequency band. At later times, Fermi acceleration be- 
comes operative when the blast Lorentz factor drops below a certain critical value, leading to 
the recovery of the standard afterglow light curve. Interestingly, the observed drop-out bears 
resemblance with the fast decay found in gamma-ray bursts early X-ray afterglows. 
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1 INTRODUCTION 

The prompt emission of gamma-ray bursts (GRB) is followed by 
an afterglow phase commonly attributed to the synchrotron emis- 
sion of shock accelerated electrons (Meszaros & Rees 1997). As 
the blast wave sweeps up matter and decelerates, the dissipated 
power decreases and the emission shifts to longer wavebands (e.g., 
Piran 2005). To model this afterglow emission, one usually en- 
codes the acceleration physics in a minimal/maximal Lorentz fac- 
tor (7miii / 7max), in the spectral index s of the electron spectrum, 
in the fraction of the dissipated energy that is carried by these 
electrons and in the fraction es stored in magnetic turbulence. 

However, our understanding of relativistic Fermi accelera- 
tion has made significant progress in the last decade, to an extent 
that motivates a direct test against observational data. The conver- 
gence of analytical calculations and extensive particle-in-cell (PIC) 
numerical calculations has led in particular to the following pic- 
ture. At (superluminal) ultra-relativistic shock waves, Fermi power- 
laws cannot develop because the particles get advected to the far 
downstream along with the magnetic field lines to which they are 
tied (Begelman & Kirk 1990), unless strong turbulence has been 
excited on scales significantly smaller than their Larmor radius 
(Lemoine et al. 2006; Niemiec et al. 2006; Pelletier et al. 2009). 
In very weakly magnetized shocks, such turbulence can be excited 
by micro-instabilities in the shock precursor and therefore Fermi 
acceleration can develop, as confirmed by recent PIC simulations 
(Sironi & Spitkovsky 201 1). The critical level of magnetization be- 
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low which this turbulence develops depends on the shock Lorentz 
factor (Lemoine & Pelletier 2010, 2011) as indeed, such instabil- 
ities can grow only if their growth timescale is shorter than the 
timescale on which the unshocked plasma crosses the shock pre- 
cursor and, the stronger the upstream background magnetization, 
or the larger the shock Lorentz factor, the shorter the precursor. 

In practice, one may expect Fermi acceleration to proceed un- 
hampered if the blast wave propagates in a weakly magnetized ex- 
ternal medium such as the interstellar medium (ISM). In magne- 
tized stellar winds, however, one might expect to see signatures of 
the above microphysics of Fermi acceleration, all the more so at 
early stages when the blast Lorentz factor is large. Such signatures 
would open a window on the physics of collisionless relativistic 
shocks as well as on the astrophysics of GRB afterglows. This mo- 
tivates the present study, which proposes to compute the afterglow 
light curve of a gamma-ray burst propagating in a magnetized stel- 
lar wind from "ab initio" principles regarding Fermi acceleration. 

The recent studies of Li & Waxman (2006) and Li (2010) of- 
fer an interesting perspective on this problem. From the observa- 
tion of X-ray afterglows on sub-day scales, these authors infer a 
strong lower bound on the upstream magnetic field of gamma-ray 
bursts afterglows, > 200 fiG r\/^ (no the upstream density in 
cm~^); Li (2010) actually derives a significantly stronger bound 
by considering on equal grounds the long lived high energy emis- 
sion > 100 MeV. This implies that either micro-instabilities have 
grown and excited the magnetic field to the above values, or the pre- 
existing magnetic field itself satisfies this bound. While the former 
is expected if the circumburst medium is ISM like, the latter corre- 
sponds to a magnetized circumburst medium. It is this possibility 
that will be addressed and tested in the present work. 
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Section [2] presents an analytical discussion of the model and 
is followed by numerical calculations of the light curve in Sec. [3] 
Section |4] discusses the results in the light of modem X-ray after- 
glows. 



2 PHYSICAL MODEL 

We consider the following fiducial values for the parameters char- 
acterizing the afterglow. The ejecta is composed of a homoge- 
neous shell of width A — cT with T = lOTi s in the station- 
ary frame, with (isotropic equivalent) bulk kinetic energy E — 
10^''i754erg and Lorentz factor 7oj — 3007cj,2.5- This outflow 
impinges on a stellar wind with density profile /9„ = Ar^^, 
with A — 5 X 10^^ At g/cm and (toroidal) magnetic field profile 
B„ = lO^B. G (r-/10^^cm)~^ the variables A, and B, en- 
code our uncertainty on the density and the magnetic field. The 
value B, — 1 corresponds to the formation of a magnetar with 
surface field ~ 10^^ G after the collapse of a WRq progenitor 
star. The magnetic field of Wolf-Rayet stars, which are considered 
as potential progenitor stars for gamma-ray bursts, are not known, 
but surface values as large as 1 — 10 x 10'^ G have been consid- 
ered (Ignace et al. 1998). For these parameters, the magnetization 
cTw = B^/ (47rpwC^) ^ 1.8 xlO"*i?2A7^ it does not depend on 
r. Note that the central assumption here is that of a relatively high 
magnetization of the external medium; the density profile does not 
play a crucial role and similar effects can be observed in a constant 
density medium of sufficient magnetization, as discussed briefly in 
Sec.il 

The proper density of the ejecta riej = 
Eej/ (47rr^7ej AmpC^) ; therefore the density contrast be- 
tween the ejecta and the external medium is not very large, 
(riej/nw)'/^ ^ SlT-'/SoiVs-Bsi^^^'^", implying that the 
reverse shock propagates at relativistic speeds in the ejecta (Sari & 
Piran 1995). The shocked material - which we denote as the blast 
- thus moves with initial Lorentz factor 

7b ^ ^ f ^ 110 T-'^'EltA:'^' . (1) 

It remains constant as long as the reverse shock is crossing the shell 
(Sari & Piran 1995, Beloborodov & Uhm 2006). Approximating 
the velocity of the reverse shock as c in the ejecta frame, the reverse 
shock has crossed the outflow at radius — ^'^cT (in the station- 
ary frame), corresponding to observer time tx — 5 sec (1 -f z)T-i, 
with z the GRB redshift. 

Beyond r-^, the blast Lorentz factor decreases according to 
7b — 7b, X {r /rx)'^^"^ in the adiabatic regime, 7b, x correspond- 
ing to Eq. lO. The relationship between observer time, radius 
and blast Lorentz factor then becomes (fobs/ix) — i^/^x)^ — 
(7b/7b,x)"''- 

We now come to the modelling of the electron population in 
the blast. Following Lemoine & Pelletier (2010, 2011), we define 
the parameter Yinsst = ^i^^o'vjJsh^ which characterizes whether in- 
stabilities may develop or not in the shock precursor, hence whether 
Fermi cycles can develop or not. The parameter ^b denotes the 
fraction of incoming matter energy through the shock 47b pwC^ 
that is carried by the accelerated and returning particles (i.e. the 
beam). By returning, it is meant those incoming protons that are re- 
flected on the shock front, which constitute an essential ingredient 
of the shock formation. These reflected protons exist even in the 
absence of Fermi powerlaws. Through mixing with the unshocked 



plasma, these returning particles (along with the accelerated parti- 
cles) induce two-stream or filamentation micro-instabilities in the 
shock precursor, on scales close to the electron to ion skin depth 
c/cjpc — c/cjpi. The filamentation instability has time to grow 
only if Yinst ^ 1, while other two stream instabilities may grow 
faster but are inhibited once the background electrons are heated 
to relativistic temperatures in the shock precursor (Lemoine & Pel- 
letier 20 11). For this reason, we consider only the growth of the fil- 
amentation instability in the following. We define a threshold value 
Yc such that if Yinst < Yc, micro-instabilities can grow and al- 
low Fermi cycles to develop, as discussed further below, while if 
Yinst > Yc, instabilities cannot grow, hence Fermi cycles do not 
develop. 

One must expect Yinst > Yc in the early stages of the after- 
glow, since 

Yinst 43 B?A,-''"Tr'/'£;5\/''Cb/_i (^-^^ ' , (2) 

with ay — for r < rx and ay ~ —1 for r > rx. PIC sim- 
ulations indicate that ^b,-i = Cb/0.1 ~ 1 (Siroiri & Spitkovsky 
2011). 

Early on, as Yinst > Yc, micro-instabilities are quenched 
by advection of the plasma through the shock front, hence the 
magnetic field is everywhere transverse to the shock normal 
without substantial inhomogeneity on short scales. In this case, 
Fermi acceleration cannot develop as particles are advected 
with the magnetic field lines to the far downstream. Never- 
theless, the electrons acquire part of the kinetic energy of the 
incoming protons in the shock transition (as viewed in the 
shock frame). A detailed understanding of this process is still 
lacking but current PIC simulations confirm the above, even 
in the absence of filamentation in the precursor. In particular, 
Sironi & Spitkovsky (2011) observe that reaches the value 
of 0.1 at a magnetization au = 10~^, for 7sii ~ 20 and larger. 
We adopt this value in the following. For simplicity, we model 
the shock heated electron distribution as a restricted power- 
law with 7 

max — - S'Yniin • 

The minimal Lorentz factor 7niin is 
then related to through 7min = ee7b('Tip/me)as with as = 

[{S - 2)/{s - 1)] [1 - (7ntax/7min)^- = ] [l - (7max/7min ) 

a normalization prefactor of order unity, which depends (slightly) 
on the modelling of the energy distribution; we adopt s — 2.4, 
an ad-hoc choice here as well motivated by simplicity (i.e. s will 
not change once Fermi acceleration becomes effective). Although 
the electrons are heated in the shock transition, the magnetic 
field is only compressed, so that the magnetic field in the blast 
frame Bb = 47b -Bw- In terms of the conventional parameter es 
describing the fraction of energy carried by the magnetic field in 
the blast, es = i3b/(87reb) = 2a„. 

As the blast Lorentz factors decreases beyond rx, so does 
Yinst, until Yinst < Yc eventually. The filamentation instabil- 
ity now has several to many e— folds of growth times before the 
plasma is advected through the shock front. This has several con- 
sequences of importance. First of all, the upstream electrons are 
heated in the micro-turbulence in the shock precursor (Spitkovsky 
2008, Lemoine & Pelletier 2011) and they therefore reach rough 
equipartition with the incoming protons after the shock transition, 
as observed in PIC simulations (Sironi & Spitkovsky 2011). This 
implies ee — 0.5. Furthermore, a micro-turbulent magnetic field 
is generated on skin depth scales up to es of a few percents. We 
adopt es = 0.05 as a fiducial value in what follows, in the absence 
of more detailed results from PIC simulations in the parameter 
range of interest. Finally, as discussed above, the micro-turbulence 
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unlocks the particles off the magnetic field lines and allow them 
to scatter repeatedly back and forth the shock wave, leading to a 
power-law extension beyond the (relativistic) thermal population. 
This fact has been clearly observed in the PIC simulations of Sironi 
& Spitkovsky (2011), for 7sh — 20 and upstream magnetization 
10~^ (and mass ratio mp/rrie — 16). Note that the same simula- 
tions at magnetization 10~* do not observe signs of Fermi acceler- 
ation, suggesting that Yc < 0.5. In the following, we keep manifest 
the dependence on Yc. To model the resulting electron distribution, 
we use a power-law between 7min and 7max, with 7,niii related to 
as before (although here, te — 0.5); we keep s = 2.4. This implies 
that we do not distinguish between the thermal and the powerlaw 
tail populations; this is a good approximation, since both radiate 
in synchrotron, hence the above simplification only affects the flux 
normalization at high energies. 

Synchrotron energy losses provide an upper 
bound on the maximal Lorentz factor, 7max 4 x 

^06^1/8^3/8^-3/8^-1/4^(^/^.^)3/4 ^. > ^-^^ 

es,-i.3 = es/0.05). In the present case, the actual limiting 
factor for 7max comes from the scattering properties of accelerated 
particles in the micro-turbulence, as discussed in Pelletier et al. 
(2009). Indeed, Fermi cycles can develop if the angular scattering 
in the micro-turbulence dominates over the regular Larmor orbits 
in the background field on a cycle timescale, which requires 
r^,o ^ £c{5B / B)^ , with 7-l,o the Larmor radius in the background 
field 47b -Bw, £c the micro-turbulence scale and SB the micro- 
turbulence strength. Both simulations (e.g., Sironi & Spitkovsky 
2011) and analytical arguments (e.g. Lemoine & Pelletier 2011) 
indicate that the relevant length scale is the ion skin depth Si (as 
measured in the upstream frame), while {5B / B)'^ ~ eB/cr^. This 
implies a maximal Lorentz factor 

7max ~ S)xW'Eli*T-^'^AV*eB,-i.zB:^ , (3) 

assuming here r > r x . 

Several remarks are in order at this stage. We do not consider 
the issue of the evolution of the micro-turbulence in the down- 
stream, which remains an open problem in this field (Gruzinov 
& Waxman 1999, Medvedev & Loeb 1999). The typical Larmor 
radius of electrons of Lorentz factor 7bmp/me ~ 7min reads 
tl ~ Si/e'g'^ , hence the first generations of accelerated electrons 
only probe the vicinity of the shock front in terms of Si, where 
the turbulence should not have evolved strongly. In our case, the 
electron population develops on a a dynamic range ~ 5B/B ~ 
(es/cTw)^^^ 10^, see Eq. (jSj. Therefore the highest energy elec- 
trons explore the blast up to ~ 10*(5i (given that the scattering 
length in the micro-turbulence scales as 7^); admittedly, one cannot 
exclude that the turbulence evolves on such length scales. For refer- 
ence, the ion skin depth 5i ~ 2.8 X 10^ craEl{^Tl''^A-^{r/ry,). 

We also neglect the influence of extra large scale sources of 
turbulence, associated with e.g. instabilities of the blast itself (e.g. 
Levinson 2010), or with the interactions of the shock with inhomo- 
geneities of the wind (e.g. Sironi & Goodman 2007). This is jus- 
tified insofar as the strong background magnetic field effectively 
prevents particles located further than ~ tl.o away from the shock 
front to return to the shock front, and tl.o ~ SijaX!'^ is already 
much smaller than the typical scales at which such instabilities de- 
velop. This means that particles that undergo Fermi cycles cannot 
experience turbulence sourced on scales larger than tl.o- Finally, 
the present study does not discuss the impact of pair loading in 
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Figure 1. X-ray light curve of a gamma-ray burst located at z = 1, with 
fiducial parameters as described at the beginning of Sec. |2] in solid line, 
result of the model, revealing the flux drop-out at tobs ~ 100 s and the 
recovery at late times ~ Iff* s. In dotted line, the same model, assuming 
however that micro-instabilities can grow at all times in front of the forward 
shock (thus implying = 0.5, = 0.05 at all times); this situation also 
corresponds to what would be seen for magnetizations (Ju < 10"®, all 
other parameters remaining unchanged. 



front of the shock wave (e.g. Beloborodov 2005, Ramirez-Ruiz et 
al. 2007), which will be addressed in a future work. 



3 LIGHT CURVE 

The above description provides the necessary ingredients to com- 
pute the light curve. We will be mostly interested in the X-ray 
afterglow, which probes the highest energy electron population 
at the early stages of the afterglow. We rely on the model intro- 
duced by Beloborodov (2005), which assumes that electrons cross 
the shock front, get instantaneously accelerated to a powerlaw, 
then cool adiabatically and through synchrotron / inverse Comp- 
ton losses. This model fits nicely the present description and the 
present hierarchy of timescales: tacc ^ ^loss ^ A7b/c and 
rL,o ^ A7b (the blast width in the blast rest frame). We have 
added the spectral contribution of fast cooling electrons to the 
model of Beloborodov (2005) in order to discuss the X-ray light 
curve. We also take into account inverse Compton cooling follow- 
ing the parametrization of Sari & Esin (2001) as discussed in Li & 
Waxman (2006): in particular, at early times when the blast magne- 
tization es = 2(Jw ^ and electrons cool rapidly, the Compton 
parameter Y\c — (ee/es)^''^, while at late times, in the slow cool- 
ing regime, Yic ~ 1. The radiative loss time is then written as 
tsyn/(l + Yic), with tsyn the synchrotron loss time in the blast 
frame. At rx and beyond, Klein-Nishina effects are not significant 
since /ii^min7miii ~ 7bmeC^ for the fiducial values (at rx) and 
7max ~ 37min before recovery. The deceleration of the blast wave 
is followed by solving the equations of the mechanical model of 
Beloborodov & Uhm (2006). The ejecta and the blast are assumed 
homogeneous and once the reverse shock has crossed the ejecta, 
its contribution is discarded from the equations of motion. We also 
assume an adiabatic evolution of the blast wave. This is clearly jus- 
tified at early times, when = 0.1; at late times, = 0.5 but the 
emission takes place mostly in the slow cooling regime, therefore 
this remains a reasonable approximation. 
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Figure 2. Upper panel: evolution in time of the blast Lorentz factor (left 
y— axis) and of the instability pai'ameter l^nst (right j/— axis); the onset of 
Fermi acceleration occurs when Yinst < 1 in the model shown in Fig. [T] 
corresponding to t^hs — 10^ s. Lower panel: evolution with time of the 
minimal and maximal observer frame electron synchrotron frequencies; at 
times < 10'* s, Fermi acceleration does not take place, a powerlaw cannot 
develop hence i-'max is 1 order of magnitude larger than the frequency 
^min- Note that i-'min E^its the X-ray band around to^s ~ 100 s. At late 
times, Fermi powerlaws develop and the range i/min — fmax broadens sig- 
nificantly, leading to the recovery of the standard light curve. 



Figure[T]presents the resulting light curve in the energy inter- 
val 0.3— 10 keV. The parameters correspond to the previous fiducial 
values; we have also adopted es = 0.05, Yc = 1 and z = 1. AX 
early times, tobs ^ 10* s, Yinst > Yc hence there is no Fermi 
power-law, only a thermal electron population extending over half 
an order of magnitude, implying that the synchrotron emission ex- 
tends over an order of magnitude. The (observer frame) frequency 
^min associated to 7min reads 

0.2 eSb 



1 + Z TTleC 
1.6 X 10^* Hz 



1 + 2 



l/2y-3/2 



(4) 



with a„ = —1 for r < r x , Oi/ = —3 otherwise. Consequently, for 
r > rx, meaning tohs > 5s(l + z)Ti, the minimum frequency 
drops rapidly out of the X-ray band. This is accompanied by a dras- 
tic reduction in flux as the maximal frequency i/max (~ 9fmin in 
the absence of Fermi powerlaw) also exits progressively out the X- 
ray domain. Given the strong dependence of Umin on r, the drop-out 
occurs shortly after tx'- in detail, defining the drop-out time td- 
as that at which Vmhi — 0.7 x 



lO" Hz, 



110 s Eli' a:'''' BT'et^Ml + z 



,1/3 



(5) 



This timescale does not depend on the duration of the prompt emis- 
sion (although it cannot of course be shorter). The shape of the 
light curve during the drop-out is affected by our assumption of 
a restricted powerlaw; a more detailed modelling of the electron 
spectral distribution (e.g. Giannios & Spitkovsky 2009) is required 
to refine the prediction for the light curve in this region. One should 
also account for the delay associated with emission away from 
the ligne of sight or from within the blast, which would lead to 
a smoothing of the light curve on a timescale ~ fobs at fobs- 

The above simplified model predicts no flux in the X-ray band 
between the completion of the drop-out, roughly a factor of a few 



beyond td-o and the recovery, i.e. the time at which Finst = Yc- 
This latter timescale free can be written 



tr 



0.9 X 10SS54A7^B^y,"^Cb,-i(l + 



(6) 



At free, Fermi cycles develop on a very short timescale compared 
to the dynamical timescale, hence emission can take place up the 
maximal frequency i^max corresponding to 7max, with 



7.6x 10'* Hz E7Ml'-^^B~^Y-'e'l^^_ 



i.asb.-i ( - — 

'•rcc 



-3/2 



(7) 

The strong dependence of this maximal frequency on the parame- 
ters suggests that a variety of effects could take place; in particu- 
lar, one might observe a weak recovery in the X-ray band or even 
no recovery at all. Caution has to be exerted however, since when 
tohs ^ 10* s the Lorentz factor of the blast has dropped to mod- 
erate values ~ 20, hence additional effects may come into play. In 
particular, one cannot rule out the emergence of new instabilities 
at scales larger than Si that would push 7max hence Vmax to much 
larger values. At even later times, jet sideways expansion affects the 
dynamical evolution; as the above one dimensional model ignores 
such effects, we stop the calculation alt = W s. 

Figure |2] summarizes the evolution with fobs of the main pa- 
rameters and allows to understand better the behaviour of the X- 
ray light curve shown in Fig. [T] Note that one does not expect a 
drop-out in the optical, since !^max ^ 10*^ Hz for the present fidu- 
cial values and Uynin crosses the optical shortly before free- Such 
a drop-out could only be seen if free were made much larger than 
10* s, e.g. by increasing the magnetization, see Eqs. (|4]l,([5},(|6j. 



4 DISCUSSION 

Using "ab initio" principles of relativistic Fermi acceleration now 
tested in extensive PIC shock simulations, we have calculated the 
X-ray afterglow light curve of a GRB propagating in a magne- 
tized stellar wind of magnetization ctu ~ 10^*, assuming oth- 
erwise standard GRB parameters. We have shown that the inhi- 
bition of relativistic Fermi acceleration in magnetized shocks at 
high Lorentz factor leaves a distinct signature in the light curve, 
in the form of a fast drop-out shortly after the end of the prompt 
emission, around fobs ~ 100 s, with a recovery at late times 
~ 10* s. The latter depends more strongly on the model pa- 
rameters, in particular magnetization so that one may envisage a 
variety of situations beyond that described: e.g., no drop-out if 
(Tu < 10^*' - ceteris paribus - or a drop-out with no recov- 
ery 

if cTu > 10"^ Although we have calculated the light curve 
for a stellar wind profile, similar effects might be observed for 
a constant density circumburst medium, provided the magnetiza- 
tion is large enough. In particular, one would observe a drop-out at 



180 sec (1. + z)*/^B 



2/3 ipl/3 -1/3 4/3 
-3 ^54, 



for a magnetic 



B-3G and density n — no cm , followed by 

^ 10* sec (1. + z)B%'Eli\^'^'Y-'/'^-'Jf. 



td- 

field B = 10" 
recovery at tree 

Therefore the present results extend beyond the stellar wind case 
and may be applicable to both long and short GRBs. 

Interestingly, recent Swift observations have revealed a rather 
complex early X-ray afterglow light curve in a subset of GRBs, 
with a fast decay at fobs ~ 100 s followed by a form of plateau that 
joins a more standard light curve at later times > 10* s (Nousek 
et al. 2006, O'Brien et al. 2006). High latitude emission is consid- 
ered as a possible explanation for the steep decay phase, although 
modelling the plateau phase with the afterglow brings in additional 
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constraints on the overall GRB model (e.g. Panaitescu 2007). The 
present scenario could account for two of these observed features - 
the initial fast decay and the late time recovery - but it does not ex- 
plain the emergence of the plateau. The following briefly addresses 
these issues in turn. 

Regarding the fast decay phase, the present scenario predicts 
an exponential decay and a clear spectral transition from hard to 
soft as the peak of the emission exits the X-ray band. One does 
not therefore expect a perfectly smooth transition from the prompt 
emission to the fast decay phase. As discussed in Sec. [3] addi- 
tional theoretical developments are required to provide a detailed 
light curve around 100 s. Nevertheless, it is of interest to note that 
Sakamoto et al. (2007) have reported evidence for an exponential 
decay component on top of a powerlaw decaying component in the 
early X-ray light curve. Furthermore Zhang et al. (2007) have ob- 
served a pronounced hard-to-soft spectral transition during the fast 
decay in two thirds of GRBs that show a fast decay (see also Yo- 
netoku et al. 2008). Their phenomenological model involves an ex- 
ponentially cut-off energy spectrum, the peak energy Ec of which 
moves out of the X-ray band during the fast decay according to 
Ec oc f^^^' with Qc — 1 — 1.5. This fits quite well the present 
picture, considering in particular that Umin oc t~y^^^ [Eq. Q], We 
also note that some short GRBs show an exponentially decaying 
light curve around 100s, well beyond the prompt emission, accom- 
panied by spectral evolution, such as GRB050724 (Campana et al. 
2006), while some show fast decay without apparent late time re- 
covery (e.g. GRB051210, GRB060801, see Nakar 2007). 

Regarding the shallow decay phase, a contribution from the 
reverse shock has been envisaged in Uhm & Beloborodov (2007) 
and Genet et al. (2007), although some tuning appear required to 
ensure a smooth transition to the recovery phase. After the present 
work was submitted, a paper by Petropoulou et al. (201 1) appeared, 
arguing that the shallow decay phase can be accounted for by the 
low energy tail of the synchrotron self Compton component. Alter- 
natively, one could try to explain the shallow decay phase with an 
inefficient contribution of the forward shock, the efficiency increas- 
ing with time and reaching its maximum at recovery of the standard 
light curve. This could be accomplished if a small fraction x^e with 
X <C 1 is stored in an accelerated electron powerlaw at that time - 
beyond the thermal component that amounts to ee(l — x) ~ ~ 
with X rising up to ~ 1 at recovery. Granot et al. (2006) and loka 
et al. (2006) have proposed a similar scenario, with varying mi- 
crophysical parameters during the shallow decay phase. Our model 
assumes a sharp transition at Yinst ~ 1 between no acceleration 
(i.e. X — 0) and fully efficient acceleration (x ~ 1), but what ac- 
tually happens at Vinst ~ C(l) when a few efolds of growth of 
the turbulence can occur is not known, since we have at our dis- 
posal the results of only two simulations at au = 10^'* and 10~^. 
Moreover, one should recall that current PIC simulations probe tiny 
timescales in regards of the GRB timescales and that these simu- 
lations do not yet converge to a stationary shock state (Keshet et 
al. 2009). One thus cannot exclude that inefficient acceleration oc- 
curs at Finst ~ but goes undetected in current simulations; 
dedicated PIC simulations on long timescales appear required to 
probe this transition region. Alternatively, if the jet is structured in 
energy and Lorentz factor per solid angle, an observer may receive 
emission from regions of different Lorentz factors than that on the 
line of sight (e.g. Panaitescu 2007); if the Lorentz factors in those 
off-axis regions are such that Yinst < Yc, one might detect low 
flux emission, corresponding to x < 1 and possibly a shallow de- 
cay phase. Yet another possibility is that of a clumpy circumburst 



medium, with clumps of various sizes, provided Tc <§i r /jh El As 
the causal region of lateral extent r /7b contains many clumps, one 
does not expect a bumpy signature in the light curve. However, one 
would collect only a fraction x < 1 of the expected X-ray flux due 
to acceleration in the fraction /c < 1 of the clumps that carry a 
magnetization such that Finst < at a given time. As the overall 
density and Lorentz factor decrease, /c increases and so does x un- 
til recovery, which corresponds to Yinst Y in the smallest scale 
clumps that carry most of the mass. In each of the above scenar- 
ios, one would expect a smooth transition in the light curve with 
no spectral evolution between the shallow decay phase and the late 
time normal decay phase, as reported by Nousek et al. (2006) and 
O'Brien et al. (2006). 

More work is certainly warranted to discuss these aspects in 
more detail and to compare the properties of the light curve to ob- 
servational data in the relevant wavelength domains. One may in 
particular expect the inverse Compton GeV emission to provide 
further constraints on the present scenario. 
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